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Abstract 



We give a prove of evidence that the original power counting by Weinberg can 
be applied to estimate the contributions of the operators contributing to the vr- 
deuteron scattering length. As a consequence, vr-deuteron observables can be used 
<*' , to extract neutron amplitudes — in case of vr-deuteron scattering this means that the 

vr-neutron scattering length can be extracted with high accuracy. This result is at 
■ variance with recent claims. We discuss the origin of this difference. 

in 

o _ 

^ . 1. In absence of neutron targets, it became common practice to use few-body 

nuclei as effective neutron targets. To extract vr-neutron (yr-n) amplitudes, vr- 
Q I deuteron (vr-d) scattering has been studied in the past. This program can be 

successful only when both the proton observables and the few-body corrections 
are known to high accuracy. As the former can be measured directly, they do 
not cause any problem. For the latter the development of chiral perturbation 
theory for few-nucleon systems promised a controlled, model independent, 
I high precision evaluation of the corresponding amplitudes. This program was 

put forward in a series of publications, e.g., for vr-d scattering (see [1] and 
references therein). 

All those analyses are based on the conjecture of Ref. [2] that the transition 
operators for reactions on nuclei with external sources can be constructed per- 
turbatively within chiral perturbation theory. The resulting operators are then 
to be convoluted with the appropriate nuclear wave functions. For this to work 
it needs to be assumed that the contribution of few-nucleon counter terms to 
the transition operators can be estimated on the basis of naive dimensional 
analysis. If we apply this recipe to nd scattering, the leading counter term 
(Fig. 1 (d)) appears at 5th order — two orders down compared to the leading 
few body correction (Fig. 1 (c)). This was recently confirmed by an explicit 
calculation of the counter term contribution assuming natural strength for the 
transition operator [4]. 
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In contrast to this it was found recently that a logarithmic scale dependence 
shows up in the leading few-body correction to 7r-d scattering (Fig. 1 (c)) that 
calls for a counter term already at this very order [3,4] (see also [5]), which 
would preclude any high accuracy extraction of 7r-n scattering parameters 
from TT-d data. This finding is based on a perturbative treatment of one-pion 
exchange. 

In contrast, in this letter, we demonstrate by an explicit numerical calculation 
that the logarithmic divergence disappears, if we treat the one-pion exchange 
non-perturbatively to obtain the wave function. This explains, why previous 
studies basically lead to identical numbers for the leading few-body correction 
although very different wave functions were employed (see discussion in Ref. 
[1]) ^ . Stated differently, we will show that the contact term that necessarily 
arises at next-to-leading order (NLO), when pions are treated perturbatively 
in the wave function, can be calculated once the pion exchange is included 
non-perturbatively in the wave function. This was already conjectured in Ref. 
[5], but not shown explicitly. 

Thus the main goal of our study is to investigate the regulator dependence of 
the leading few-body correction. Since we are going to employ wave functions 
that contain non-perturbative pion contributions, this study can only be per- 
formed numerically. We will use deuteron wave functions that were constructed 
for cut-offs that vary over a wide range (A = 2-20 fm~^ = 400-4000 MeV). 
The procedure of their construction is described in Ref. [8] and will be briefly 
reviewed below. Already in Ref. [1] a mild cut-off dependence was reported 
for calculations using wave functions with non-perturbative pions, when the 
regulator was changed from 500 to 600 MeV. This might either be because of 
the absence of the logarithmic divergence due to the wave functions used or 
simply because the coefficient in front of the logarithm is accidentally small. 
Due to the large range of variation of cut-off values used here we are in the 
position to answer this question: we will show that there is no sign of a log- 
arithmic regulator dependence of the results as soon as the complete wave 
functions are used. The consequences of this observation will be discussed in 
the final section. 

In Ref. [9,10], it was stressed that care has to be taken when calculating pion 
reactions on nuclei. There it was shown that a subtle cancellation pattern exits 
between contributions from loops in one-body and few-body operators. This 
has the effect that the static pion exchange is an excellent approximation to the 
exact result for the leading few-body corrections to 7r-d scattering. Therefore 
here we will focus on the static exchange only. 

2. In our investigation we use the wave functions constructed as outlined in 

^ Please note that in Ref. [6], it was shown that in the deuteron channel pions 
should not be treated perturbatively. 
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(a) (b) (c) (d) 

Fig. 1. Typical contributions to ird scattering. Diagram (a) and (6) show the tree 
level and the one-loop contribution to the one-body term, the pion rescattering 

contribution is depicted by (c) and diagram (d) shows a two-nucleon contact term. 
In this figure solid (dashed) lines denote nucleons (pious) and ellipses the deuteron 
wave function. 

Ref . [8] . They emerge as a solution of the Schrodinger equation 

nip) = G{p) J d'p'v{p,p')h{p,p')nip') , (1) 

where G{p) — (— e —p^/M)~^ denotes the two-nucleon propagator with e and 
M for the deuteron binding energy and the nucleon mass, respectively. The 
leading order potential V{p^p') comprises contributions from both the one- 
pion exchange as well as a contact term as depicted in Fig. 2 (see [7]). As 
regulator function we use 

/a(p,P ) = exp I 1 . (2) 

For a given value of the regulator A the only free parameter is C — the strength 
of the contact term as depicted in Fig. 2(b). For this study, this parameter was 
adjusted such that the deuteron binding energy was reproduced, to exclude any 
dependence of the results on an incorrect asymptotic behavior of the deuteron 
wave function. We checked that the description of the phase shifts in the "^Si- 
^Di channel is comparable to the one obtained in [8]. This numerical study can 
only be conclusive, when wc cover a wide range of cut-offs. We decided to use 
values of A between 2 and 20 fm~^ (400-4000 MeV). This range starts below 
the chiral symmetry breaking scale of A^^ ^ 1000 — 1200 MeV and extends to 
values larger by a factor of 4. In this range, we also observe the appearance of 
spurious bound states in the ^Si-^Di channel. However, their energies are large 
and, therefore, these bound states should not affect any low energy physics. 

For comparison we also prepared a series of wave functions from only contact 
NN interactions (thus omitting diagram 2(a) in the potential). These wave 
functions are denoted by '^{p) in what follows. Again we impose a regulator 
as given in Eq. (2) in the Schrodinger equation. 

For completeness, we summarize the binding energy results and some wave 
function properties for both series of wave functions in Table 1 and 2. 
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A 


Eq 


T 


Pd 


As 


V 


rd 


Qd 


2 


2.225 


28.91 


5.24 


0.839 


0.030 


1.889 


0.3005 


4 


2.225 


45.48 


8.23 


0.866 


0.027 


1.933 


0.2827 


6 


2.224 


62.33 


6.94 


0.866 


0.025 


1.932 


0.2704 


8 


2.225 


75.95 


6.76 


0.864 


0.026 


1.926 


0.2676 


12 


2.227 


85.80 


7.14 


0.864 


0.026 


1.925 


0.2675 


16.5 


2.214 


102.50 


7.08 


0.862 


0.026 


1.929 


0.2676 


20 


2.210 


115.07 


7.07 


0.861 


0.026 


1.929 


0.2675 


Expt. 


2.225 






0.8846 


0.0256 


1.9671 


0.2859 



Table 1 

Summary of some deuteron properties obtained from ^'J^ — the wave functions with 
non-perturbative one-pion exchange for various cutoffs. Here, the cut-off A is given 
in fm~^, the binding energy and kinetic energy Eq and T in MeV, the asymptotic 
S'-state normalization As is in fm^^/^, the point nucleon radius in fm, and the 
quadrupol moment in fm^. 77 is the ratio of the asymptotic S- and D-state normal- 
ization. 
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Eq 


T 


As 


rd 


2 


2.225 


32.60 


0.76 


1.728 


4 


2.225 


69.59 


0.72 


1.617 


6 


2.225 


106.76 


0.71 


1.585 


8 


2.225 


143.99 


0.70 


1.569 


12 


2.225 


218.51 


0.69 


1.555 


16 


2.225 


293.04 


0.69 


1.547 


20 


2.225 


367.59 


0.69 


1.543 


Expt. 


2.225 




0.8846 


1.9671 



Table 2 

Summary of deuteron properties obtained from the wave functions ^ — where 
only a contact interaction was used in the potential — for various cutoffs. The nota- 
tion is the same as in Table 1. 

3. Let us now turn to the calculation of the leading few-body correction to the 
Trd scattering length. As stated in the introduction we will exclusively focus 
on the static contribution. The corresponding expression reads 

^(static) ^-^ I Sp£q^>l{p- g)^4*A® ' (3) 

J q 

where ^ = m^/(327r'^/^ (1 + m7r/(2Mjv))). Clearly, no physical quantity can 
be regulator dependent. Naive dimensional analysis does not require a two- 
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Fig. 2. Contributions to the NN potential at leading order: the one pion exchange 
(a) and a momentum independent contact term (b). 
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Fig. 3. Results for the A dependence of the leading few-body correction from Eq. 
(3). The 'x' symbols show the results of the numerical evaluation of the integral 
using the full wave functions whereas the symbols shows the result for 
wave functions with only point interactions ^, downscaled by a factor of 4. The 
dashed line is a fit to the latter of the form A + Bln{A), as described in the text. 
The solid line shows the fit of a constant to the former for values of the cut-off 
larger than 6 fm~^. 

body counter term in the same order as this first three-body correction. Thus, 
studying the A dependence of the given integral tells, whether a counter term 
of naively higher orders is needed in conjuction with this few-body correction. 
We evaluated numerically the integral in Eq. (3) using both, the wave func- 
tions from the full LO NN potential {k, = tt) and those from only the point 
interaction {k — no tt), as described in the previous section. 



The results are shown in Fig. 3. Here the x symbols emerged from the cal- 
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culation with the wave functions — where the non-perturbative pion ex- 
change was included — whereas the plus symbols stem from the calculation 
using — that does not have any pion exchange in the wave function. For 
the latter, our results clearly show the ln(A) behavior in accordance with the 
findings of Refs. [3,4,5]. To show that exphcitly, we fitted a logarithm to our 
results, which is also shown in the figure. The perfect agreement shows that 
we can recover the previous results numerically in our cut-off range. However, 
the calculation using the full wave functions shows almost no regulator de- 
pendence at all for cut-offs above ^ . Thus, as soon as the pion exchange is 
included non-perturbatively in the wave functions, no counter term is needed 
at the order of the leading few-body corrections to absorb the regulator de- 
pendence of the pion exchange contribution. For a numerical comparison with 
previous work, we fitted a constant to our results for A > 6 fm~^. In this way, 
we obtain for the static three-body contribution a**"**'^ = 2.12 m~^, which is 
in good agreement with the previous calculations [1] . 

We also checked that there is no unnatural enhancement of the leading nNN 
ttNN counter term due to the wave function at small distances: the contri- 
bution of this term to the scattering length — from an explicit evaluation of 
diagram (d) in Fig. 1 — was in line with the counting, when the transition op- 
erator was assumed to be of natural strength [11]. The same observation was 
also made in Ref. [4]. This shows that as soon as non-perturbative pions are 
included in the construction of the wave function, the wave function at the 
origin assumes natural values. Combining this finding with the observation 
made above that there is no regulator dependence of the leading few-body 
correction, there is no reason to change the original Weinberg counting. 

We conjecture that the full wave function is driven not only by the the binding 
momentum of the deuteron, as ^'^^ but also by a second scale, probably 
'iTT'TT or f^^. The numerical calculation does not easily allow to identify such 
a scale in the wave functions. For illustration, we compare the two types of 
wave functions in r-space in Fig. 4. For ^J^, we see that the wave functions 
converge for radii larger than R = 0.7 fm. For smaller r, the wave functions 
are not unique in our range of cutoffs and show different numbers of nodes 
coming from the spurious bound states. They remain, however, small for these 
distances. This is different for the point-hke wave functions, which, for r = 
increase linearly with the cut-off. We conclude that the finite range of the 
pion-exchange regularizes the wave functions in the sense that they remain 
small for small distances independent of the cut-off and oscillate. Both effects 
reduce the contribution of the small distance behavior to the scattering length. 
A more systematic insight would obviously be highly welcome. 



^ Note that we observe a mild regulator dependence of the integral for regulators 
below A-^, in line with the findings reported in Ref. [1]. 



6 



S.Oi 




r[fm] 



1.5| 




r [fm] 



Fig. 4. The upper panel shows and the lower one the s-wave of for various 

values of A. Note the different scales of the figures. 

4. To summarize, as soon as a finite range potential is used in the construction 
of the deuteron wave functions, the matrix element for the leading few-body 

correction to the vr-d scattering length gets independent of the regulator. In 
particular, the logarithmic divergence that emerges in case of wave functions 
from a point-like NN potential disappears. We have demonstrated this in a 
numerical study using wave functions constructed for a large variety of cut- 
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offs. Our numerical results extend to cut-offs, which are much larger than the 
chiral symmetry breaking scale. For the same range of cut-offs, we numerically 
observe the expected logarithmic divergence for the case of ip"^" . Therefore, 
we are confident that the range of cut-offs is large enough to take our result 
as a prove of evidence that no counter term is needed at the same order where 
the leading few-body correction appears. This is the precondition to allow for 
a high accuracy extraction of the vr-n scattering length from 7r-d data. 

Please note that the present discussion has large implications also for the 
investigations of other reactions. Based on our findings the calculations for 
7d ^ Ti^d [12,13], TT^He-^ vr^He [14], -n^d -inn [15], and 7^ ^ vr+nn [10] 
indeed were performed with the accuracy as given in the publications. 

Finally, we note that the observation, that the one-pion exchange, which leads 
to a finite range interaction, changes the divergence structure of the theory, 
is not unique to vr-d scattering. The probably most famous other example in 
this context is the three-nucleon system in the spin 1/2 channel: as long as 
only point-like two body interactions are included, a three body counter term 
needs to be promoted to leading order [16]. However, as soon as there is a 
non-perturbative pion-exchange included, this is no longer necessary [8]. 
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